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Abstract 

A theoretical model of nuclear spin conversion in molecules controlled by an 
external electromagnetic radiation resonant to rotational transition has been 
developed. It has been shown that one can produce an enrichment of spin 
isomers and influence their conversion rates in two ways, through coherences 
and through level population change induced by radiation. Influence of con- 
version is ranged from significant speed up to almost complete inhibition of 
the process by proper choice of frequency and intensity of the external field. 
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I. INTRODUCTION 



It is well known that many symmetrical molecules exist in Nature only in the form 
of nuclear spin isomers Spin isomers are important for fundamental science and have 
various applications. They can serve as spin labels, influence chemical reactions or 
tremendously enhance NMR signals PIH. Progress in the spin isomers study depends heavily 
on available methods for isomer enrichment. Although this field has been significantly 
advanced recently (see the review in 0), one needs more efficient enrichment methods. 

Recently two enrichment methods based on optically induced change of molecular level 
populations have been proposed |7]|| . The purpose of the present paper is to investigate the 
isomer enrichment caused by optically induced coherences in the molecule. 

II. QUANTUM RELAXATION 

Nuclear spin isomers of molecules were discovered in the late 1920s. The most known 
example is the isomers of H2. These isomers have different total spin of the two hydrogen 
nuclei, 1 = 1 for ortho molecules and 1 = for para molecules. Symmetrical polyatomic 
molecules have nuclear spin isomers too. For example, CH 3 F can be in ortho, or in para state 
depending on the total spin of three hydrogen nuclei equal 3/2, or 1/2, respectively (see, 
e.g., PJ). Different spin isomers are often distinguished also by their rotational quantum 
numbers. Consequently, all rotational states of the molecule are separated into nondegen- 
erate subspaces of different spin states. Schematically, these subspaces for the case of two 
spin states, ortho and para, are presented in Fig. |1| (For a moment, an optical excitation 
shown in the ortho space has to be omitted.) 

Nuclear spin conversion can be produced by collisions with magnetic particles. This is 
a well-known mechanism for the conversion of hydrogen imbedded in paramagnetic oxygen 
If a molecule is surrounded by nonmagnetic particles, collisions alone cannot change the 
molecular spin state. In this case the spin conversion is governed by quantum relaxation 
which can be qualitatively described as follows. Let us split the molecular Hamiltonian into 
two parts, 

H = H + hV, (1) 
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where Hq is the main term which has the ortho and para states as the eigen states (states 
in Fig. [l|); V is a small intramolecular perturbation able to mix the ortho and para states. 
Suppose that the test molecule was placed initially in the ortho subspace. Due to colli- 
sions the molecule starts undergo fast rotational relaxation inside the ortho subspace. This 
running up and down along the ortho ladder proceeds until the molecule reaches the ortho 
state m which is mixed with the para state k by the intramolecular perturbation V. Then, 
during the free flight just after this collision, the perturbation V mixes the para state k with 
the ortho state m. Consequently, the next collision can move the molecule to other para 
states and thus localize it inside the para subspace. Such mechanism of spin conversion was 
proposed in the theoretical paper |L0| (see also |TT| ). 

Relevance of the described mechanism to actual spin conversion in molecules is not at 
all obvious. The problem is that the intramolecular perturbations, V, able to mix ortho and 
para states are very weak. They have the order of 10 — 100 kHz (hyperfine interactions) 
which should be compared with other much stronger interactions in molecules, or with gas 
collisions. Nevertheless, the experimental and theoretical proves have been obtained that 
spin conversion in molecules is indeed governed by quantum relaxation [§,12-14]. Although, 



only three molecules (CH 3 F, H 2 CO and C2H4) have been studied in this context so far, it 
is very probable that spin conversion in other polyatomic molecules of similar complexity is 
governed by quantum relaxation too. It is useful for the following to give a few examples of 
spin conversion rates. Most studied is the spin conversion in 13 CH3F which has the rate, 

7l3 /P = (12.2 ± 0.6) • 10" 3 s _1 /1brr, (2) 

in case of pure CH 3 F gas ||. Spin conversion in another isotope modification, 12 CH 3 F, is by 
almost two orders of magnitude slower ||. Similar slow conversion was observed in ethylene, 
13 CCH 4 , JTJ, 

j eth /P = (5.2 ± 0.8) • 10~ 4 s-VTorr. (3) 

These data show that spin conversion by quantum relaxation is on 9 — 11 orders of magnitude 
slower than the rotational relaxation, v ~ 10 6 — 10 7 s^ 1 /Torr. 
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III. OPTICAL EXCITATION 



In this section we will analyze the spin conversion by quantum relaxation in the presence 
of a resonant electromagnetic radiation. The level scheme is shown in Fig. |l|. In order 
to reveal the main features of the phenomenon we will consider the process in a simplest 
arrangement. First of all, we assume that molecular states are nondegenerate and that 
only one ortho-para level pair, m — k, is mixed by the intramolecular perturbation, V. 
Monochromatic radiation is chosen to be in resonance with the rotational transition m — n. 
In this arrangement the molecular Hamiltonian reads, 

H = H + hG + hV. (4) 

New term, G, describes the molecular interaction with the radiation, 

G = —(E d/h) cosu^t, (5) 

where E and ujl are the amplitude and frequency of the electromagnetic wave; d is the 
operator of the molecular electric dipole moment. We have neglected molecular motion in 
the operator G because homogeneous linewidth of pure rotational transition is usually larger 
than the Doppler width. 

Kinetic equation for the density matrix, p, in the representation of the eigen states of 
the operator Hq has standard form, 

dp/dt = S-i[G + V,p], (6) 

where S is the collision integral. 

Molecules in states (m, n, k) interacting with the perturbations G and V constitute only 
small fraction of the total concentration of the test molecules. Thus, one can neglect collisions 
between molecules in these states in comparison with collision with molecules in other states. 
The latter molecules remain almost at equilibrium. Consequently, the collision integral, S, 
depends linearly on the density matrix for the disturbed states m, n, and k even in one 
component gas. Further, we will assume model of strong collisions for the collision integral. 
The off-diagonal elements of S are 

S Jj' = ~ T Pjj': j,f e m,n, fc; j ^ f. (7) 
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Here, j and j' indicate rotational states of the molecule. The decoherence rates, T, were 
taken equal for all off-diagonal elements of S. 

Collisions cannot alter molecular spin state in our model. It implies that the diagonal 
elements of S have to be determined separately for ortho molecules, 

Sjj = -vPo(j) + vw (j)p ; Po = J2Poti)'i j G ortho, (8) 

j 

and for para molecules, 

s ij = -vPpti) + v w p(3)Pp'i P P = ^2p p U)> J g P ara - ( 9 ) 

3 

Here p a , p p and w Q (j), w p (j) are the total concentrations and Boltzmann distributions of 
ortho and para molecules. The rotational relaxation rate, u, was taken equal for ortho and 
para molecules, because kinetic properties of different spin species are almost identical. 
One can obtain from Eq. (0) an equation of change of the total concentration in each 



spin space [p|| . For example, for ortho molecules one has, 

dp /dt = 2Re(ip mk V km ). (10) 

In fact, this result is valid for any model of collision integral as long as collisions do not 
change the total concentration of molecules in each subspace, i.e., J2j Sjj — if j G ortho, 
or j G para. 

Let us approximate the density matrix by the sum of zero and first order terms over 
perturbation V, 

p = pW+pW. (11) 

In zero order perturbation theory the ortho and para subspaces are independent, i.e., p^ k = 
0. Consequently, equation of change ( |T0D is reduced to, 

dp /dt = 2Re(ip% k V km ). (12) 

Note that the spin conversion appears in the second order of V. 

Kinetic equations for zero and first order terms of the density matrix are given as, 

dp ( ®/dt = S ( ® -i[G,p m ]- (13) 



dp {l) /dt = S {1) - i[G, - i[V, p (0) ]. (14) 
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We start with zero order perturbation theory. In this approximation, the para subspace 
remain at equilibrium, 

pf{j)=pfw p {j)- j G para. (15) 

Here p^ is the total concentration of para molecules. 

Density matrix for ortho molecules can be determined from the two equations which 
follows from Eqs. ©,®,(0): 

dp^(j)/dt = -vpf>{j) + vw {j)pf> + 2Re(ip^(m\n)G nm )[8 jm - 8 jn }; 
dpf\m\n)/dt = -r p V>(m\n) - iG mn [pf\n) - pf\m)\. (16) 

We will assume further the rotational wave approximation, 

Gmn = —Ge int ; G = Eod mn /2h; £1 = tui — io mn , (17) 

where the line over symbol indicate a time-independent factor. Rabi frequency, G, is assumed 
to be real. Rotational relaxation, u, and decoherence, T, are on many orders of magnitude 
faster than the spin conversion. It allows to assume stationary regime for ortho molecules, 
thus having, dp^'(j)/dt = 0. The substitution, p^>{m\n) = p^(m\n) exp (— iQt), transforms 
Eqs. ([16]) to algebraic equations which can easily be solved. Thus one has, 



Pi 0) (j)=pi 0) 
p^{m\n) =p<®iGAw 



2T G 2 Aw 



r| = r 2 + AYG 2 /u; Aw = w (n) -w (m). (18) 

We turn now to the calculation of the first order term, p^ k , which has to be substituted 
into the equation of change (|l2l). The density matrix element, p^ k , can be found from the 
two equations which are derived from Eqs. fl7D,(|14]), 

dp^Jdt + Tp™ + iG m J2 = -iV m M\k) - pfH]; 

dp^/dt + Yp { 2 + zG nm p« = iV mk pf{n\m). (19) 

Substitution, 

V , — Vp iuJt (u) = u) ,V — -n^p iuJt - nP) - rtWp'^-Wfa)* (on) 
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transforms Eq. (19) to algebraic equations from which one finds p^ k . Then, Eq. flT2|) gives 
the following equation for the total concentration of ortho molecules, 

dpo = 2lVl2R ^ + ^ + ^MKk)-p^(m)]- l Gp^\n\m)_ 



dt 1 1 F(Q) 

F(Q) = (T + ioo)[T + + u)} + G 2 . (21) 



After an appropriate change of notations, the right hand side of Eq. (|21|), coincides 
formally with the solution [|15| for the work of weak optical field in the presence of strong 
optical field. Strong field splits the upper state m on two, which appears as two roots of the 
denominator of Eq. ( pl|) being the second order polynomial on u. It results in two ortho-para 
level pairs mixed by the perturbation V instead of one pair in the absence of an external 
field. In analogy with the optical case, one can distinguish the isomer conversion caused by 
population effects (terms proportional to the level populations Pf'(k) and p^(m)) and by 
coherences (term proportional to the off-diagonal density matrix element, ^\n\m)). 

Eq. fl2l|) describes time dependence of the concentration of ortho molecules, p Q , in the 
second order of V. One can neglect at this approximation small difference between p Q and 

. The density of para molecules can be expressed through the density of ortho molecules 
as, p p ^ = n — p(°), where n is the total concentration of the test molecules. Using these points 
and zero order solution given by Eqs. ([15]), ([18]), one can obtain final equation of change for 
ortho molecules, 

dp / dt = nj p - p j; 7 = 7 op + 7 po + 7„ + -f coh , (22) 
where partial conversion rates were introduced, 

P -I- -4- 1, A 

7op = 2\V\ 2 w p (k)f(n); /(fi) = Re- 



Fin) 

lpo = 2\V\ 2 w (m)f(n); 



,- l2 G 2 Aw „ T-ifi 

7 - = 2m fuW Re TW (23) 

Here the rates 7 op , 7 po , and j n are due to molecular level populations and the rate 7^ is 
due to coherences. The terms 7 op and 7 po are the only ones which remain in the absence of 
an external field. In the radiation free case (G = 0) the result (p3|) becomes identical with 



the solution given in [11 



IV. ENRICHMENT 



Solution to Eq. (p2] ) can be presented as, p Q = p + 5p exp(— jt), where time-independent 
part is given by 

- Po = n^. (24) 
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Without an external radiation (at the instant t — 0), the equilibrium concentration of para 
molecules is equal to, 

w n (m) n 

^--"•' o »° \ (t )tU 'i < 25 > 



For simplicity, the Boltzmann factors in Eq. fl25|) were assumed to be equal, w p (k) = 
w (m) = which implies, 7 op = 7 po . External field produces a stationary enrichment 
of para molecules. One can derive from Eqs. (123), (E3), 



(3 = J^L- -1 = 1-2^. (26) 

An enrichment coefficient, (3, is defined here in such a way that (3 = if there is no external 
electromagnetic field. Enrichment of ortho molecules is equal to —(3. Note, that the en- 
richment, (3, does not depend on the magnitude of intramolecular perturbation V. It is the 
consequence of the assumption that only one ortho-para level pair is mixed. Enrichment, (3, 
depends on the ratio of Boltzmann factors, w (n)/w (m), but does not depend on the mag- 
nitude of w (n) itself. In further numerical examples relative difference of the Boltzmann 
factors will be chosen as w (n) = 1.2w. 

One needs to specify a few other parameters in order to investigate properties of the 
optically induced enrichment. We will use, where it is possible, parameters relevant to the 
spin conversion in 13 CH3F. Thus the decoherence rate, T, will be chosen equal 6 MHz, which 
corresponds to the gas pressure of pure CH 3 F equal 0.2 Torr |J. Rotational relaxation, z/, 
will be chosen by one order of magnitude slower than the decoherence rate, v = 0.1T. 

Expressions for the enrichment, (3, are given in the Appendix. If ortho-para mixing is 
performed for a degenerate pair of states m — k (u = 0), the enrichment of para states, (3, 
has one peak at Q = 0. This peak is determined mainly by population effects. 

More interesting is the case of nondegenerate states (u 7^ 0). In this case one has two 
peaks, at Q = —uj and at Q = (Fig. ||). Peak 1 (Q = —u) is due to the coherent effects 
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determined by ■y co h- Peak 2 (Q = 0) is mainly due to the optically induced level population 
changes determined by j n . In the case of well-separated peaks (T B <^ u), the amplitudes of 
the peak 1 and peak 2 read, 



A 1 
A 2 



Aw G 2 
w T 2 + G 2 



+ 



D 



D = T 2 + 5G 2 /4 - TG 2 /u + u 2 /4; 



Aw G 2 
"f 1 



B 



r 3 



D 



T 2 + G 2 + cu 2 



(27) 



Amplitude of the peak 2 grows rapidly with G up to j3 ~ 4.5%. Amplitude of the peak 1 
grows with G to even bigger value (3 ~ 5.5%. These data are shown in Fig. |3| (upper panel) 
where points are obtained by fitting an exact expression (|26D by two Lorentzians and solid 
curves are given by Eqs. (|27|). u was chosen equal 130 MHz which corresponds to the 
ortho-para level gap in 13 CH 3 F 0. Note, that there is no optically induced enrichment if 
Aw = 0. 

Widths of the enrichment peaks are given by the expressions, 



Wt = 2VT 2 + G 2 ] W 2 = 2T 



B- 



f28) 



The two enrichment peaks experience completely different power broadening which are 
shown in Fig. |3] (low panel). Peak 1 has the width much smaller than the Peak 2. Solid 
curves in Fig. |3| (low panel) are given by Eqs. ([28|). Points are obtained by fitting the exact 
expression for enrichment, (3, by two Lorentzians. 



V. CONVERSION 



Conversion rate in the presence of an external electromagnetic field has complicated 
dependence on radiation frequency detuning, Q, and Rabi frequency, G. It is convenient to 
characterize the conversion rate in relative units, 

Irel = —-l, (29) 
T/ree 

where 7/ ree is the field free conversion rate. Similar to the enrichment, (3, this parameter does 
not depend on magnitude of the perturbation V and on absolute values of the Boltzmann 
factors. 

Expressions for the conversion rate, 7 re ;, are given in the Appendix. In the case of 
degenerate ortho-para level pair m — k (u — 0), 7 re ; has narrow negative structure at Q = 
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(Fig. UJ). Amplitude of this dip grows rapidly with increasing G. If G ^> T ^> u, and all 
Boltzmann factors have the same order of magnitude, the conversion rate at Q = is given 
by 

Irel ~ £2 " 1- (3°) 

If G is large, the relative conversion rate, 7 re j ~ —1, which corresponds to 7 ~ 0. Thus the 
spin conversion can be inhibited by radiation having large G and Q = 0. 

In the nondegenerate case (u 7^ 0) conversion rate, j re i, has two peaks (Fig. [5]). If these 
peaks are well resolved, Tb ^ uj', conversion rate is given by, 

lrel ~ v 2 + (n + uo'f + w r| + ^ U 2) ' { } 

where new parameters are determined as, 



r'-r(i + f i^); ^4-^]. (32) 



The two peaks in the conversion have Lorentzian profiles and amplitudes determined by 
the expressions, 

G 2 . Aw G 2 fT 1 



These amplitudes are shown in Fig. |6| (upper panel) by solid curves. Points are obtained by 
fitting the exact solution by two Lorentzians. In this examples, u was chosen equal 130 MHz. 
Peak 4 at Q = is proportional to the ratio of Boltzmann factors. This peak does not grow 
significantly with G. The peak 3 at Q = —uj 1 almost does not depend on Boltzmann factors 
and at large G grows up to 7 re z = (u/T) 2 . Thus strong electromagnetic field can speed up 
the conversion significantly, viz., by two orders of magnitude in our numerical example. 
Widths of the two peaks in the conversion rate are determined by the equations, 

W 3 = 2V; W A = 2T B} (34) 

and have very different field dependences. In fact, the peak 4 at Q = is broadened even 
faster than 2Tb- Peak 3 at Q = —uj has almost no power broadening, if G uj (Fig. |^, low 
panel). Solid curves in Fig. Bl low panel, corresponds to the Eqs. (0). Points are obtained 
by fitting the exact solution by two Lorentzians. 
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VI. CONCLUSIONS 



We have shown that an external resonant radiation can influence spin isomer conversion 
in two ways, through level populations and through optically induced coherences. The coher- 
ences introduce in the process new and interesting features. In many cases, the coherences 
play more important role than the level populations. This analysis have been performed in 
a simplest arrangement in order to reveal the main features of the phenomenon. 

Optically induced coherences introduce extra resonances both in enrichment and in con- 
version frequency dependences. These new resonances are important for future experimental 
realizations of the optical control of isomer conversion. First, they give convenient oppor- 
tunity to find coincidences between molecular transitions and available sources of powerful 
radiation. Second, an observation of the phenomenon will be easier also because electro- 
magnetic radiation can significantly speed up the conversion. Thus steady state enrichment 
can be achieved much faster than without field. It allows to work at low gas pressures where 
described above effects can be achieved at smaller radiation intensity. 

Another advantage is that one can use resonance at which there is no large radiation 
absorption and, consequently, no significant level population change by radiation. It should 
decrease some spurious effects, like molecular resonance exchange. The latter effect can 
cause serious problems in realization of the optically induced enrichment by population 



Here we give expressions for the enrichment, (3, and conversion rate, ^ re \. Enrichment of 
para molecules is given by Eq. ( p6|) . For equal Boltzmann factors, w (m) = w p {k) = w, one 
has 



effects 0| . 



VII. APPENDIX 




(35) 



and an approximate expression in case of small enrichment, 



7n + Icoh 



(36) 



Using Eqs. (^), this expression can be reduced to, 
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p 



Aw 



1/4 



Y/u - 3/4 



+ 



D 



y 2 + g 2 + (q + u) 2 Y 2 B + n 2 (r| + fi 2 )[r 2 + G 2 + (fi + cu) 2 ] 

D = Y 2 + 5G 2 /4 - YG 2 jv + uj 2 /A. 



(37) 



There are two peaks in enrichment, at O = — uj and at O = 0. These peaks have Lorentzian 
shape if uj 3> Yb- In the degenerate case (u = 0) there is one peak of complicated form at 

n = o. 

The spin conversion rate, 7 re /, is defined as, 



7reZ 
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- 1. 



7/r 



(38) 



Using Eqs. 



7reZ = 9 



one can obtain after straightforward calculations the following expression, 

[l + g)(uj 2 -Y 2 ) + 2uj{Sl + uj') 

Y' 2 + (n + uj') 2 ^ 

Aw g 2 /r i 



w r| + Q 2 V i/ 2 



l + g 



'l + g)(uj 2 -Y 2 ) + 2u(n + uj') 

r 2 + (n + uj') 2 



+ 



Aw G 2 (l+g)(2T 2 + uj 2 )-uj( y Q + uj') 
2w f| + fi 2 F 2 + (fi + cj') 2 

G 2 



i" = r(i + 



uj' = uj(1 — 



(39) 



Y 2 + uj 2 ' 

Thus the conversion rate has two peaks situated at Q = —uj' and at — 0. Position of the 
peak at Q = uj' depends on the field intensity. 

In the case of degenerate ortho-para level pair {uj = 0), frequency dependence of 7 re ; 
has complicated shape with a dip in the center at Q = (Fig. @). The whole structure is 
described by the expression, 



irel 



[1 + g)G 2 + Aw G 2 



Y' 2 + n 2 



w Y 2 B + tt 2 



[l+g)G 2 \ + Y'Y 



Y' 2 + n 2 



Y' 2 + n 2 



(40) 
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FIG. 1. Molecular ortho and para states. Bent lines indicate the rotational relaxation. Vertical 
line shows an optical excitation. V m k gives the ortho-para state mixing by an intramolecular 
perturbation. 
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FIG. 2. Enrichment of para molecules, (3, as a function of radiation frequency detuning, $7, 
= 130 MHz and two Rabi frequencies, G = 1 MHz and G = 6 MHz. 
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FIG. 3. Amplitudes of the peaks in enrichment, (5 (upper panel) and width (FWHM) of these 
peaks (low panel) for the case of w = 130 MHz. 
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FIG. 4. Relative conversion rate, 7 re j, for the case of degenerate ortho-para states, uj = and 
the two values of Rabi frequency, G = 1 MHz and G = 1.5 MHz. 
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IG. 5. Relative conversion rate, j re i, for the case of nondegenerate ortho-para states 
130 MHz) at two Rabi frequencies, G = 1 MHz and G = 6 MHz. 
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FIG. 6. Amplitudes of the conversion pesks, ^reU 

(upper panel) and widths (FWHM) of these 
peaks (low panel) for the nondegenerate case, oj = 130 MHz. 
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